Deterministic implementation of weak quantum cubic nonlinearity 
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We propose a deterministic implementation of weak cubic nonlinearity, which is a basic building 
block of a full scale CV quantum computation. Our proposal relies on preparation of a specific 
ancillary state and transferring its nonlinear properties onto the desired target by means of deter- 
ministic Gaussian operations and feed-forward. We show that, despite the imperfections arising 
from the deterministic nature of the operation, the weak quantum nonlinearity can be implemented 
and verified with the current level of technology. 

PACS numbers: 



Ever since it has been first mentioned by Feynman [l| , 
quantum computation has been the holy grail of quan- 
tum information theory, because the exponential speedup 
it offers promises a significantly better way to tackle cer- 
tain non-polynomial computational tasks. However, in 
order to achieve this, high-order non-linear resources are 
needed. The original approach to quantum computing 
relied on manipulation of quantum systems with a lim- 
ited Hilbert space 0], but one can expand this concept 
to arbitrary quantum systems and say that the universal 
quantum computation is defined by its ability to emu- 
late a Hamiltonian in form of an arbitrary (Hermitian) 
polynomial of annihilation and creation operators Q . 

It is a great advantage of the continuous variables (CV) 
quantum systems [3] that Hamiltonians composed of first 
(linear) and second (quadratic) powers of quadratic op- 
erators x and p can be readily implemented. However, it 
can be easily checked that these are not enough and that 
an access to cubic Hamiltonian with at least third power 
of quadrature operators is necessary. Unfortunately, the 
currently achievable experimental interaction strengths 
are too low compared to noise to be of use. 

Fortunately, the need for currently unavailable cubic 
unitary evolution may not be so dire. Let us recall the 
original statement of Lloyd and Braunstein If one 
has access to Hamiltonians A and B, one can approxi- 
matively implement operation with Hamiltonian i[A, B]. 
Approximatively is the key term here, meaning that the 
desired operation is engineered only as a quadratic poly- 
nomial of the interaction time: 
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Consequently, even the initial operations need not to be 
unitary - their quadratic approximations are fully suffi- 
cient. What this means is that if we take interest in a 
sample cubic interaction with Hamiltonian H cx x 3 , we 
need not to implement the unitary e lxx , but it is enough 
to be able to perform operation 



This is the lowest order expansion for which the commu- 
tator trick (p} works, but let us start with the real lowest 
order expansion, 1 + /3x 3 , where j3 is a complex number. 
This expansion behaves as weak cubic coupling if beta 
is imaginary and has the added benefit that is can be 
used to compose ^ when the respective values of (3 are 
complex and chosen properly. In principle, even this gate 
can be further decomposed into series of 1 + jx (7 6 C) 
operations @ . These phase sensitive gates can be imple- 
mented probabilistically on a traveling beam of light by 
subsequent application of photon subtraction and pho- 
ton addition, represented by operators a = (x + ip) / \f2 



and |g-|g]. They are very useful for preparing various 
ancillary states, but for use in a full-fledged information 
processing we are interested in their deterministic imple- 
mentation. 

Ideal implementation: To this end we employ the 
approach of [9| , thoroughly discussed in [lj| , where it was 
suggested that unitary operation acting on a state can 
be deterministically implemented with help of a proper 
resource state, a QND coupling, a suitable measurement 
and a feed- forward loop. Explicitly, for operation 0(x) 
acting on pure state = J ip{x)\x)dx the resource state 
is 0{x)\p — 0). After QND coupling, represented by 
unitary £/qnd(A) = e 1 ^ 2 ^ 1 , is employed and the overall 
state is transformed to 



ip(x)0(y)\y - Xx, x)dxdy, 



(3) 



the ancillary resource mode gets measured by a homo- 
dyne detection. We can for now assume A = 1, as the 
overall message remains unchanged. For any detected 
value q the output state is 



i/j(x)G(x + q)\x)dx. 



(4) 
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To obtain the desired result, one either post-select only 
for situations when q = was detected, or apply a feed- 
forward which would compensate for x + q in the argu- 
ment of the operator. It has been shown in Q that if 



2 



the desired operation 0(x) is a unitary operation driven 
by a Hamiltonian of order n, the feed-forward operation 
requires Hamiltonian of order n — 1. Explicitly, imper- 
fections in operator <D(x + q) = exp[ix(x + q) 3 } can be 
compensated by unitary operator [/ff = exp[— ix(3qx 2 + 
3q 2 st)}, which is a combination of displacement, squeezing 
and phase-shifts. The operation @ we are interested in 
is not unitary, but since it is an approximation of unitary 
driven by a cubic Hamiltonian, a feed-forward of squeez- 
ing and displacements should perform adequately, up to 
some error. We'll get to this issue later. In fact, the oper- 
ations available for feed-forward limit us in what we can 
do. With squeezing and displacement we can implement 
only cubic operations. Of course, with them we could 
also tackle Hamiltonians of the fourth order, and so on. 
And there is another limitation - since the feed-forward 
must be deterministic and noiseless, and therefore uni- 
tary, it can be only used to deterministically compen- 
sate unitary (at least approximatively) operations whose 
Hamiltonian is hermitian. Therefore we cannot use the 
trick with implementing a series of 1 + jx operations, we 
have to implement operation ([2]) in one go. Consequently, 
we need a sufficiently complex resource state. 

Resource state generation: Let us now shift our 
attention to the required resource state. In realistic, 
even if idealized, considerations one has to, instead 
of a position eigenstate, use a squeezed state S\0) = 
[j exp(— x 2 /g)\x)dx}/(ng) 1 / 4 , which approaches the ideal 
form as g — > oo . The resource state can now be expressed 
as: O(x)S\0) = SO(x/y/g)\0) which is a state finite in 
Fock basis with superficial layer of squeezing. As it has a 
finite structure, the state can be engineered by a sequence 
of six photon additions [ll[ or photon subtractions [Hj]. 
This is an extremely challenging task, let us therefore 
first focus at the lowest nontrivial cubic hamiltonian ex- 
pansion, Os(x) = 1 + xx 3 , which is a feasible extension of 



recent experimental works [Tj 
state looks as 



S(l + X 'x 3 )\0) =S 



The appropriate resource 



(5) 




with x' 
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This state can be generated from a 



squeezed state by a proper sequence of photon subtrac- 
tions and displacements [l3| . which acts as (a — a) (a — 
(3)(a — 7)5*10). Since the squeezing operation transforms 
the annihilation operator as S^aS — fxa — i>a\ where 
H = cosh(ln v /g) and v = sinh(ln»/g), the required dis- 
placements can be obtained as a solution of set of equa- 
tions: 



A = afry, a + 13 + 7 = 0, 
2V2V 3 = A X ', iv 2 + 3/w = (a/3 + a 7 + h) > (6) 

where A is a constant parameter related to normalization. 



The solution exists and it can be found analytically as 



7 = -(a + /3). 



(7) 



Here £ and £ are solutions of the set of equations 

xy + C x = 0, y - x 2 - C 2 = 0, (8) 



where d = v 3 2V%X , C 2 = iv 2 + 3/xza The solu- 
tions of © always exist and they can be obtained ana- 
lytically using the Cardan formula. The squeezing used 
in the state generation can be in general different from 
the squeezing in ([5]). However, squeezing can be con- 
sidered to be a well accessible operation and we shall 
therefore not deal with this in detail. It should be noted 
that an alternative way of preparing the state ([5]) lies in 
performing a suitable projection onto a single mode of 
a two-mode squeezed vacuum state. Engineering of the 
proper measurement, which too requires three APDs and 
three displacements, leads to similar equations as in the 
previous case (JB]) with solution of the same form. 

Realistic implementation: With the resource state 
at our disposal we can now look more closely at the two 
ways to implement the gate, the probabilistic and the 
deterministic, in order to compare them and see, what is 
the manifestation of high order nonlinearity in the deter- 
ministic case. The probabilistic implementation is rather 
straightforward. Using the resource state ([5]) we are able 
to transform the initial state to 
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ip{x)0(x)e' 



-x 2 /2g 



x)dx, 



(9) 



and as the squeezing of the resource state approaches 
infinity, the produced state approaches its ideal form. 
The final state is always pure and the actual compo- 
sition of the operator O n {x) can be arbitrary, allowing 
us, for example, to implement operator O in n different 
non-unitary steps. On the other hand, if the resource 
squeezing is insufficient compared to the distribution of 
the state in phase space, it seriously affects some proper- 
ties of the state - for example, moments of x quadrature 
may not be preserved any more. 

But let us move towards the more interesting part, 
the deterministic approach. In this case the operation 
produces a mixed state 



P = \ P{q)\%)mdq. 



(10) 



Here, P(q) represents the probability of measuring a spe- 
cific outcome q and 
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^x)e^ x+q ^/ 29 0{x + q) e -*x9 3 -*3x(*?W 9 ) \ x } dXj 



(11) 
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where Mr is the norm of the resource state, stands for 
the respective quantum state corrected by feed-forward. 
Ideally, 0(x + g) e - i3 x( x « 2+x2 9) O(x), but this relation 
can obviously work only when both x and q are small 
enough for the exponent to be reasonably approximated 
by the finite expansion O n . It is therefore quite unfortu- 
nate that the very condition required for the operation to 
work flawlessly, the need for g — > oo, is compatible with 
the feed-forward only in the limit of x 0. To quantify 
these properties in greater detail we need to employ a 
suitable figure of merit. 

Analysis: To evaluate the quality of the approximate 
operation is not a straightforward task. If wc want to 
conclusively distinguish the cubic type nonlinear interac- 
tion from a Gaussian one, we can take advantage of the 
known way the quadrature operators transform: x — > x, 
p — > p + xx 1 . If wc apply the operation, in form of a 
black box, to a set of known states, we can analyze the 
transformed states to see whether the operation could be 
implemented by a suitable Gaussian, or if it is more of 
what wc aim for. The analysis can be as easy as checking 
the first two moments of the quadrature operators, be- 
cause the nonlinear dependance of (p) on (x) can not be 
obtained by a Gaussian operation, unless we consider a 
rather elaborate detection-and-feed-forward setup, which 
would, however, introduce an extra noise detectable ei- 
ther by checking purity of the state, or by analyzing 
higher moments (x 2 ) and (p 2 ). 

The case with purity of one is straightforward to verify 
- as soon as the first moments have the desired form, 
(£') = (x) and (p') = (p)+x(x ), we can be certain a form 
of the desired non-Gaussianity is at play. In the presence 
of noise, the confirmation process is more involved and 
we shall deal with it in a greater detail. 

It needs to be shown that, in comparison to the deter- 
ministic approximation, no Gaussian operation can pro- 
vide the same values of moments (x'}, (x' 2 ), and (p') 
without also resulting in a significantly larger value of 
moment (p' 2 ) caused by the extra noise. The complete 
Gaussian scheme consists of an arbitrary Gaussian inter- 
action of the target system with a set of ancillary modes 
followed by a set of measurements of these modes yield- 
ing values which are used in suitable feed-forward to fi- 
nalize the operation. In the case where the approximate 
transformations approach the ideal scenario, i.e. when 
(£') = (x), (x' 2 ) = (x 2 ), and (p') = {p)+ X {x 2 ) 1 only a 
single ancillary mode is sufficient, the optimal Gaussian 
interaction is in the QND interaction with a parameter 
A, and after a value of £ is measured by homodyne detec- 
tion, feed-forward displacement of k£ 2 ensures the correct 
form of the three moments. In the end, the Gaussian ap- 
proximated state can be expressed as 

Ps = J dxx (12) 

D s (kx 2 ) a (x\U qnd (\)ps <8> |0)(0|[>t ND (A)|x)^( K x 2 ), 

where the subscripts S and A denote the signal and the 
ancillary mode, respectively. The high order classical 
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FIG. 1: (Color online) (a) First moments relative to real part 
of a. Solid red and blue lines represent the ideal values of (p) 
and {£}, respectively. Red and blue crosses than show these 
values for deterministic non-Gaussian approximation, while 
red and blue circles do so for the Gaussian approximation. 
Dashed green line is a quadratic fit for (p) . The experimental 
parameters are g — 1 and % — 0.03. (b) Second moments 
relative to real part of a. Solid red and blue lines represent 
the ideal values of (p 2 ) and (x 2 ), respectively. Red and blue 
crosses than show these values for deterministic non-Gaussian 
approximation, while red and blue circles do so for the Gaus- 
sian approximation. 



nonlincarity is induced by the nonlinear feed-forward, 
represented by the displacement D(a). The strength of 
the QND interaction A remains a free parameter over 
which can the procedure be optimized to obtain the 
best approximation characterized by the minimal pos- 
sible value of the extra noise term in (p' 2 ). 

We analyze the aforementioned properties over a set 
of small coherent states a with |a| < 2. We compare 
the state obtained by the approximative cubic interaction 
with the state created by the Gaussian method. In prin- 
ciple, this could be done for both the deterministic and 
the probabilistic approach, but since the probabilistic ap- 
proach has the potential to work perfectly, we shall keep 
to deterministic methods in our comparative endeavors. 
For each coherent state and its cubic-gate transformed 
counterpart we can, from knowledge of the first moments 
of quadrature operators, estimate the actual cubic non- 
linear parameter and use it to construct the benchmark 
Gaussian-like state (|T2"j) . The final step is to compare the 
extra noise present in p quadrature - if the added noise for 
the approximate state is below the Gaussian benchmark 
we can assume a non-Gaussian nature of the operation. 

As an example, let us look at a particular scenario, 
in which the deterministic cubic gate was applied to a 
set of coherent states with imaginary part of the com- 
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FIG. 2: (Color online) Schematic experimental setup of the 
deterministic x 3 gate. BS-beam splitters, PBS-polarization 
beam splitters, PS-phase shift, PR-polarization rotator, 
APD-avalanche photo diode, HD-homodyne detection, D- 
displacement. 



plex amplitude constant. The effect of the operation is 
illustrated in Fig. [TJ Fig. QJi shows the first moments 
and reveals that for this purpose, effective cubic nonlin- 
earity of XcS = 0.1 can be reliably obtained for both 
the non-Gaussian and the Gaussian approach. Differ- 
ences arise, though, for the second moments, where the 
value of Gaussian quadrature moment (p 2 ) is observably 
higher than the value of its non-Gaussian counterpart. 
The values of the Gaussian moment were obtained by 
optimizations of (|T2j) for each particular value of Rea, it 
is therefore a stronger benchmark than a universal Gaus- 
sian operation, working over the whole range of Rea, 
would be. And it is still beaten by the imperfect deter- 
ministic non-Gaussian method with no squeezing in the 
ancillary mode. 

Experimental setup (Fig. [2]) The resource state 
generation requires three-photon subtraction from a 
squeezed, or a two-mode squeezed light, with appropri- 
ate displacements. Two photon subtractions have been 
already implemented [l2| and three of them are within 
reach. The resource state is then coupled with the in- 
put using a QND gate with offline squeezing [3, Il5j |. 
which can be modified as to reliably manipulate the non- 
Gaussian resource state [ItJ ■ The final step lies in per- 
forming a sequence of feed-forwards driven by homodync 
measurement of the ancilla Q. Of those, the only non- 
trivial one is given by unitary e , where the actual 



value of A depends on the measurement. This opera- 
tion can be decomposed into a sequence of phase shift 
by 1; squeezing with gain gf and another phase shift 
by 02, where the parameters satisfy: tan0i tan 02 = — 1, 
tan 0i = .9/, and (1 — gj) cos0i sin 02 = 2g/A. Adjusting 
the squeezing gain on the fly can be done by exploit- 
ing the universal squeezer [pi . [rH ] , where the amount of 
squeezing is controlled by changing the ratio of the beam 
splitter, which can be done by a sequence of polariza- 
tion beam splitter, polarization rotator, and another po- 
larization beam splitter, where the rotator controls the 
splitting ratio. The nonlinear dependance of the feed- 
forward parameters on the measurement results requires 
a sufficiently fast data processing, but that too is avail- 
able today [H] . 

Conclusion: We have proposed an experimentally 
feasible way of deterministically achieving weak nonlin- 
carity of the third order. The procedure effectively engi- 
neers the operation on a single photon level and then de- 
terministically cut'n pastes the properties onto the target 
state. However, the nonlinearity is still composed only 
of a limited number of photon, so it can be faithfully 
applied only to target states which are sufficiently weak. 
Furthermore, since there is no such thing as free lunch, 
subsequent use of the transformed state in attempts to 
generate higher nonlinearities as per Q requires higher 
and higher numbers of single photons used in the engi- 
neering. 

The approach is not flawless. There are several sources 
of noise which can be simultaneously reduced only in the 
limit of infinitely small (read unobservable) interaction. 
This is due to the finite photon approximation of the 
cubic gate not being unitary and therefore not perfectly 
correctable by the unitary feed-forward. Nevertheless, we 
have shown that even with this noise a demonstration of 
decisively no n- Gaussian high order quantum determin- 
istic nonlinearity going well beyond classical attempts, 
based on higher order nonlinearity in the feed-forward 
loop, can be observed already now. 
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